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I. INTRODUCTION 



o 

<N 

(3JQ[ Position-dependent- mass (PDM) particles have received interests over the last few decades From the very 

^ . beginning, parallel to the applicability of the PDM-settings in different areas of physics such as the many-body 

<^ 1 problem semiconductors [3[, quantum dots quantum liquids He- Clusters jg], etc. a mathematical difficulty 
has grown which is known as the ordering ambiguity. Von-Roos has, firstly, reported the problem by proposing the 

■ revolutionary form of the PDM kinetic energy operator which reads as Q 



T = — (m a pmP pm 1 + rrPprnPpm^ . (1) 



s 

Here m is the position dependent mass of the physical particle and a, (3 and 7 are some real parameters which satisfy 
r*j ■ the constraint a + /? + 7 = — 1. In the constant mass setting, (1) reduces to the standard kinetic energy operator 
d ' which is given by 
£3 " 

We note that in Eq. (1), although a + /3 + 7 = — 1 but the exact values of a, and 7 are not known and without any 
harm to the mathematics of the problem, this von-Roos parameters may accept any real number. The nature of the 

■ problem turns to the fact that the momentum operator p does not commute with the position operator x in quantum 

' theory. 

In the literature, there exist several suggestions for the von-Roos ordering parameters, for instance, the Gora's 
and Williams' ( (3 = 7 = a = —1) Ben Daniel's and Duke's (a = 7 = 0, j3 = —1) @, Zhu's and Kroemer's 
( a = 7 = -1/2, £ = 0) |l0||, Li's and Kuhn's (a = 0, /3 = 7 = -1/2) [ll|, and the very recent Mustafa's and 
£N ■ Mazharimousavi's ( a = 7 = —1/4, (3 = —1/2) [HI- It has been observed that the physical admissibility of a 
given ambiguity parameters set very well depends not only on the continuity conditions at the abrupt heterojunction 
^ ■ boundaries but also on the position-dependent-mass form. The general consensus is that there is no unique and 
universal choice for these ambiguity parameters. 



II. CONTINUES ORDERING PARAMETERS 



The idea which we shall expand in the sequel, is to construct a quantum kinetic operator which is a superposition 
of all possible values for the von-Roos parameters a, /3 and 7. In doing so, we also introduce a weight function which 
gives the distribution of the different ordering parameters. The idea can be developed in two different directions: i) 
a continues distribution and ii) a discrete distribution. Here in this section we study the continues distribution and 
in the next section, the discrete distribution will be given. 

Following to our proposal we introduce the kinetic energy operator T for a position dependent mass particle as 

f = — J J dAp {a, (3) (m a pm l3 pm' r + m? 'pm^ 'pm a ) (3) 
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in which p (a, f3) : R 2 — > R is the weight function, is the surface element on a — (3 plane and 

A = ffdAp(a,P) (4) 

can be called a normalization constant (see Fig. 1 and 2). 

We note that the third parameter 7, is not free and therefore, the integral is taken only over two of the parameters 
(i.e., a and (3). Furthermore, one may impose additional conditions on these parameters to specify the domain of the 
ordering parameters, for instance, a, /?, 7 € [— 1,0] . Latter constraint looks to be reasonable because, the well-known 
ordering-parameter sets available in the literature, belong to this domain. In spite of this fact, we keep the parameters 
free, to get any values in R, but under a rather flexible constraint. 

Let's consider the domain of a, f3 and 7 the same and fully symmetric which we called them IZi, 1Z2-.. (Fig. 2). In 
this figure, IZi is located on the plane a + j3 + 7 = —1 in the space of a, (3 and 7. In addition to TZi we also introduce 
1Z to be the projection of TZi on the plane of a — /?. Finally, we consider IZi, an equilateral triangle which leads to 

TZ = \ (5) 
l-(l + § + a)</?<§ 

in which — | < 6 E R and it shows the size of the region TZi . 
By considering the natural unit, i.e., h— 1, one finds 

m a pwPpm? + m 7 pm^pm a = (6) 

-2((l + a + f3 + a/3 + a")^- 1 -(l + ^^-^d x + -dl), 
\ y ' m A 2 m z m z m J 



which yields 



f = -i ((1 + „ + , + + » (1 + £ - £«. + i*) . (7) 



A. A UNIFORM DISTRIBUTION 

As an example, we study the case of p(a,(3) = 1 which is a uniform and symmetric distribution. By setting 
p(a,f3) = 1 in Eq. (7) it leads to 

(»> + 4, + 28) '4 - IQ + i (2*4 - i*) . ,8) 
6 V 16 v ' m A m z 2 V m z ml 



It is remarkable to observe that for V6 £ [— |, 00) non of the first two terms are zero. We note that among the possible 
values for b, b = — | corresponds to 

a = /3 = 7 = -i. (9) 

This set of parameters consequently implies 

~ -l/5m' 2 lm"\ 1 /m' 1 „ 2 \ 

Another interesting case is when 6 = 0, which makes a, /3 and 7 lie in the interval [—1, 0] . The kinetic energy after 
setting 6 = reads 

~ -1 f7m' 2 m"\ 1 /to' 1 o2 \ . , 

= "6" \ 4 TO^ _ TO 2 " J 2 ( TO 2 " _ TO J ' 



B. A NON-UNIFORM DISTRIBUTION 

Next we give an example with the non- uniform distribution function p(a,/3) which is given by 

""•«=(^TII + 0FTI)- (12) 

This is a symmetric distribution with respect to a and j3. Consequently the form of kinetic energy becomes 

1 / to' 2 to" m' 1 „,\ 

T =-o [Vi— +V2— ~— d x + -d 2 x ), (13) 

where 

(36 2 + 126 + 20) In [ 4 ( 2 +y 2 ) ) + 2 (6 3 + 36 2 + 126 + 4) (tan" 1 § + tan" 1 (1 + 6)) + 156 2 + 406 + 20 

, (14) 

4(2+2h+b 2 ) 



and 



48 (1 + 6) (tan- 1 § + tan" 1 (1 + 6)) - 24 In 



(6 + 4) In ( 4(2 t+ b f A + 2 (36 + 2) (tan" 1 § + tan" 1 (1 + 6)) - (36 + 2) 

Vo = 7 x (15) 

16 (1 + 6) (tan' 1 \ + tan' 1 (1 + 6)) - 8 In ( ^g^ ) 



One can check that 



and 
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lim r} 2 = ~ (17) 

M-§) + 3 

which are exactly as we expected for a = = 7 = — ^ . Also for the case of 6 = one gets 

7T + 10 (1 - In 2) 7r + 2-41n2 

Vi = —7T< — „, , — - 0.589 66, ri 2 = -—r —— = -0.33746. (18) 

11 6(7r-21n2) ' 12 4(7r-21n2) v ; 

III. DISCRETE ORDERING PARAMETERS 

As we mentioned before, there are some well-known ordering-parameters sets [||-[l2j which have been used widely in 
the literatures. These different sets which have been introduced for various physical problems, could be an indication of 
having a discrete ordering parameters. In this line let's introduce the kinetic energy operator with discreet distribution 
as 

X) da,P {jn a pm^ pm 1 + m? 'fimP 'pm a ^ 
f= a -^- _— , (19) 

4 2^ ^a,P 
a,P 

in which C a ^ is a real number which can be called distribution number. We again recall that 7 = — 1 — a — /3 and in 
the summation 7 is not a free index. Very similar to the case of continues distribution one may rewrite 

[(l + a + l3 + a(3 + a 2 )^-^(l + (3)^-^d x + ^d 2 x 
?=-\- ■ ( 2 °) 
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To see how this may work, we consider an equal possibility for all known ordering given before, i.e., i) Gora and 
Williams ( /? = 7 = 0, a = —1), ii) Ben Danial and Duke (a = 7 = 0, /3 = —1), iii) Zhu and Kroemer (0 = 7 = —1/2, 
= 0), iv) Li and Kuhn (a = 0, /3 = 7 = -1/2) and v) Mustafa and Mazharimousavi (a = -1/4, /3 = — 1/2 = -1/4). 
Considering these values into (20) and setting C a ^ = 1, we find 



(21) 



and consequently 

1 /43 m' 2 3m"\ 1 /m' 1 



T ~ 10 I 16 2m^) + 2[r^ dx rrfi ' ' (22) 



IV. CONCLUSION 



In this Letter, we have considered the long standing ordering ambiguity problem associated to the von-Roos kinetic 
energy operator for PDM particles. A superposition of the possible values of the von-Roos parameters has been 
used to construct a general kinetic energy operator. We have also introduced the distribution function which can be 
continues or discreet. We believe that this kinetic energy operator can assist the theoretical physicists to adjust their 
results along the line of the relevant experiments. In the other words there would be more free parameters in the 
theoretical calculations. 
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Figure Captions: 

Fig. 1: The plane oi a + (3 + j = — lin the space of a/3-y. The surface element on this plane is given too. 
Fig. 2: The possible symmetric regions of a, f3 and 7 on the plane of a + (3 + 7 = —1. The projection of these 
regions on a/3— plane is give by (5). 
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